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ABSTRACT 

We consider a newly-born millisecond magnetar, focusing on its interaction with the dense stellar 
plasma in which it is initially embedded. We argue that the confining pressure and inertia of the 
surrounding plasma acts to collimate the magnetar's Poynting-flux-dominatcd outflow into tightly 
beamed jets and increases its magnetic luminosity. We propose this process as an essential ingre- 
dient in the magnetar model for gamma-ray burst and asymmetric supernova central engines. We 
introduce the "pulsar-in-a-cavity" as an important model problem representing a magnetized rotating 
neutron star inside a collapsing star. We describe its essential properties and derive simple esti- 
mates for the evolution of the magnetic field and the resulting spin-down power. We find that the 
infalling stellar mantle confines the magnetosphere, enabling a gradual build-up of the toroidal mag- 
netic field due to continuous twisting. The growing magnetic pressure eventually becomes dominant, 
resulting in a magnetically-driven explosion. The initial phase of the explosion is quasi-isotropic, 
potentially exposing a sufficient amount of material to 56 Ni-producing temperatures to result in a 
bright supernova. However, if significant expansion of the star occurs prior to the explosion, then very 
little 56 Ni is produced and no supernova is expected. In either case, hoop stress subsequently colli- 
mates the magnetically-dominated outflow, leading to the formation of a magnetic tower. After the 
star explodes, the decrease in bounding pressure causes the magnetic outflow to become less beamed. 
However, episodes of late fallback can reform the beamed outflow, which may be responsible for late 
X-ray flares. 

Subject headings: gamma rays: bursts — magnetic fields — pulsars: general — stars: magnetic fields 
— stars: neutron — supernovae: general 



1. INTRODUCTION 

In the classical collapsar scenario for long-duration 
gamma-ray bursts (GRBs), the core of a rotating massive 
star collapses to form a black hole, whereas the overly- 
ing stellar material possesses enough angular momentum 
to form an accretion disk that persists for at least sev- 
eral seconds, long enough for its jet to breakout from 
the star (Woosley 1993; Paczynski 1998; MacFadyen & 
Woosley 1999). This accretion disk — black hole system 
then acts as a central engine for the GRB. The power 
for the explosion comes both from accretion energy, re- 
leased via neutrinos and perhaps via a magnetic mech- 
anism (e.g., the magnetic tower mechanism as proposed 
by Uzdensky & MacFadyen 2006), and from the black- 
hole rotational energy, released via the Blandford-Znajek 
(1977) mechanism. 

In the present paper we investigate an alternative sce- 
nario, in which the central object formed as a result 
the of core collapse is not a black hole, but rather 
a rapidly-rotating (millisecond) magnetar with a large- 
scale poloidal magnetic field of the order of 10 15 G. Such 
a strong magnetic field can be produced, for example, 
by a turbulent a — £1 dynamo driven by convection in 
a proto-neutron star (PNS) subject to neutrino cooling 
(Duncan & Thompson 1992; Thompson & Duncan 1993). 
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An alternative possibility is that the progenitor core of 
about 10 4 km has a magnetic field of order 10 9 G, similar 
to the field levels actually observed in some white dwarf 
of similar size. When such a highly-magnetized core col- 
lapses into a neutron star of 10 km radius, flux freezing 
leads to amplification of the magnetic field to 10 15 G, 
as discussed in Uzdensky & MacFadyen (2006). In addi- 
tion, calculations by Akiyama et al. (2003) have shown 
that turbulent dynamo driven by the magneto-rotational 
instability (MRI) in the collapsing differentially-rotating 
core is capable of producing 10 15 — 10 16 G fields within 
about 100 km from the center, on the timescale of just a 
few tens of milliseconds after the bounce (see also Ardel- 
jan et al. 2005). Whatever the origin of the very strong 
magnetic field in the PNS is, in this paper we shall take 
it for granted. Our main goal will be to investigate the 
role such a strong field plays in the explosion dynamics. 

The idea of using a millisecond magnetar as a central 
engine for gamma-ray bursts has been first proposed by 
Usov (1992) in the context of accretion-induced collapse 
of a highly-magnetic (10 9 G) white dwarf and, indepen- 
dently, by Duncan & Thompson (1992). It has been fur- 
ther developed and applied to different explosion scenar- 
ios by several authors (e.g., Thompson 1994; Yi & Black- 
man 1998; Nakamura 1998; Spruit 1999; Wheeler et al. 
2000, 2002; Ruderman et al. 2000; Lyutikov & Bland- 
ford 2002; Thompson et al. 2004; Metzger et al. 2007). 
The present paper is also devoted to investigating the 
millisecond-magnetar scenario, but viewed within the 
overall context of a collapsing star. 

At the most basic level, the main idea is that a GRB 
(or a supernova) explosion is powered by the magnetic 
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extraction of rotational energy of the newly-born rapidly- 
rotating magnetar. This magnetic luminosity operates 
alongside the much stronger neutrino cooling, which is 
the main avenue for releasing the gravitational binding 
energy of the young, still-contracting neutron star. How- 
ever, most of the neutrinos escape to infinity without 
sharing their energy with the stellar-envelope gas (un- 
less their spectrum is strongly modified by coronal pro- 
cesses, see Ramirez-Ruiz & Socrates 2005). Magnetic 
fields, on the other hand, couple to the gas tightly and 
this makes them a very efficient explosion agent. Ener- 
getically, magnetic GRB models are usually quite plau- 
sible. For example, assuming a typical surface magnetic 
field -B* = 10 15 G, a rotation rate O* = 10 4 sec -1 , and a 
radius i?* = 10 km, it is easy to see that the resulting ba- 
sic energetics and timescales fall just in the right ballpark 
to make the millisecond magnetar a plausible candidate 
for a GRB central engine (e.g., Thompson 1994). In- 
deed, the total rotational energy of a millisecond-period 
neutron star is E rot ~ 5 • 10 52 erg, which is more than 
enough to drive a long-duration GRB. The time scale 
for the energy extraction can be estimated by divid- 
ing this available energy by the total magnetic luminos- 
ity (the spin-down power), L ma g n . The latter can be 
roughly estimated by the usual pulsar luminosity for- 
mula Lmagn ~ -E>* Rt ^* c_3 : which for the above param- 
eters yields ~ 3 • 10 50 erg sec -1 , corresponding to the 
characteristic timescale of order 100 sec. Thus, from the 
point of view of the overall energetics and timescales, the 
millisecond-magnetar central engine is just a scaled-up 
version of the Ostriker & Gunn (1971) model for regu- 
lar pulsar-powered supernovae (with the magnetic field 
scaled up by three orders of magnitude and the time 
scaled down by six orders of magnitude). 

It has to be noted, however, that the plausible overall 
energetics and timescales are, by themselves, not suffi- 
cient for making a good GRB central engine model. This 
is because there are some extra physical requirements 
mandated by observations. In particular, to make a suc- 
cessful GRB, the central engine has to be capable of pro- 
ducing an energetic outflow that is (1) ultra-relativistic; 
(2) highly-collimated; and (3) baryon-free. The pioneer- 
ing works cited above have focused mostly on the ener- 
getics and timescales, but not on the mechanisms for pro- 
ducing an outflow that satisfies these requirements (see, 
however, Wheeler et al. 2000 and Bucciantini et al. 2006 
for a discussion of collimation). Also, most of these previ- 
ous models, with the notable exceptions of Wheeler et al. 
(2000) and Arons (2003), see § 12.31 have considered a 
magnetar in isolation; that is, they have completely ig- 
nored the effect of any surrounding stellar gas on shaping 
the outflow. This may be a good approximation for the 
accretion-induced collapse of a white dwarf, but it is not 
appropriate in the collapsing-star scenario. 

In contrast, in this paper we stress that the infalling 
stellar gas is still present during the explosion and needs 
to be taken into account. Thus, an important new ele- 
ment that distinguishes our model from those previous 
works is the consideration of the interaction between a 
newly-born magnetar and the stellar plasma in which 
it is initially embedded. Specifically, we argue that the 
pressure and inertia (i.e., the ram pressure) of the sur- 
rounding stellar gas acts as a natural collimator forcing 



the magnetized outflow into two tightly beamed jets. It 
also plays a crucial role in magnetic extraction of rota- 
tional energy from the magnetar. 

In order to illustrate these ideas we introduce 
the "Pulsar-in-a-Cavity" problem as a basic-physics 
paradigm for this scenario. We describe this problem 
in detail in § [5J We first give a general description of 
the problem and its various versions. Then, in § 12. li 
we consider the simplest special case of a rotating force- 
free magnetosphere inside a fixed rigid cavity. In that 
section, we first demonstrate that differential rotation of 
the magnetic field lines is inevitably established inside 
the cavity, even if the pulsar itself is rotating uniformly; 
as a result, a strong toroidal magnetic field gradually 
builds up. We then study the long-term evolution of the 
magnetic field inside the cavity and show that the mag- 
netic luminosity increases with time. We also show that 
a massive, non- force-free plasma strip unavoidably arises 
in the equatorial plane beyond the light cylinder. In § 12.21 
we discuss the subtle issue of hoop-stress collimation and 
argue that external confinement and differential rotation 
are two important ingredients for collimating relativistic 
Poynting-flux dominated outflows. We then consider, in 
§ 12.31 the case of a magnetosphere surrounded by a cavity 
with a fixed external pressure (instead of a fixed radius). 

In § |3] we discuss a specific example relevant to the 
core collapse of a massive star: a cavity formed behind 
the stalled bounce-shock at the center of the collapsing 
star. The radius of the shock stays roughly stationary on 
the timescale for magnetic fields in the cavity to grow. At 
the same time, both the ram pressure of the gas falling 
onto the cavity and the neutrino energy deposition inside 
it decrease with time. We therefore argue that at some 
point, a fraction of a second after bounce, the magnetic 
field will inevitably start to dominate the force balance, 
leading to a magnetically-driven explosion. 

In § 21 we further explore some of the astrophysically- 
interesting aspects of our model. Thus, in § 14.11 we dis- 
cuss the possibility of MHD instabilities (e.g., kink) de- 
veloping in the twisted magnetosphere and their impli- 
cations for our model. In § 14.21 we address an important 
issue of 56 Ni production and argue that the two-phase 
nature of the explosion in our model is well-suited to ex- 
plain a large amount of 56 Ni inferred from observations. 
In § 14. 31 we briefly discuss the possibility of restarting the 
GRB engine by the fall-back of the post-explosion ma- 
terial. In § 14.41 we describe an extension to our model: 
a "magnetar-in-a-tube" , motivated by the fact that the 
material along the rotation axis does not experience a 
centrifugal barrier and hence falls onto the PNS faster. 
In section § 14.51 we discuss the implications of our model 
for pulsar kicks. Finally, in § 14.61 we suggest some di- 
rections for future numerical simulations of this problem. 
We draw our conclusions in § [5l 

2. THE PULSAR-IN-A-CAVITY PROBLEM 

In order to understand how a millisecond magnetar 
central engine operates in the collapsar context, it is 
first necessary to consider the following basic physics 
problem: what happens when an axisymmetric pulsar 
is placed inside a conducting cavity filled with a low- 
density infinitely-conducting plasma (see Fig. [I])? Specif- 
ically, we are interested in a situation where the cavity 
radius Rq is much larger than the pulsar light-cylinder 
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Fig. 1. — Aligned pulsar inside an infinitely-conducting spherical 
cavity of radius Ro at t = 0. The vertical dashed lines represent 
the pulsar's light cylinder of radius -Rlc < ^0- 



radius Rlc- We call this idealized problem the Pulsar- 
in-a-Cavity problem (Uzdensky & MacFadyen 2006) and 
we propose it as the first essential step in building up a 
physical understanding of the problem. It is a modifica- 
tion of the famous problem of an axisymmetric rotating 
magnetic dipole in free space, considered by Goldreich 
& Julian (1969) as a model for an isolated pulsar's mag- 
netosphere. In some aspects, it is also similar to the 
system considered by Kardashev (1970), Ostriker (1970; 
unpublished), and by Ostriker & Gunn (1971). 

We would like to point out that adding a cavity makes 
the problem, in a sense, less fundamental since the be- 
havior of the system depends, in general, on the assumed 
physical properties of the cavity. To make the situation 
less arbitrary, we shall fix the electromagnetic proper- 
ties of the cavity by assuming that its walls are perfectly 
conducting, as is the plasma that fills the cavity. We 
shall also assume that all the field lines close back to the 
pulsar inside the cavity, i.e., that there are no field lines 
connecting the pulsar to the cavity wall. This choice 
is most natural when the bulk of the magnetic flux has 
been produced by a dynamo operating inside the neu- 
tron star and then emerged through its surface (as is in 
the case of magnetars), as opposed to a situation where 
the pulsar field lines connect directly to the wall (i.e., to 
the outer stellar envelope, as considered, e.g., by Karda- 
shev 1970 and by Goldreich et al. 1971). At the same 
time, we are still left with a lot of freedom regarding the 
mechanical properties of the cavity. Thus, we are dealing 
not with one unique problem, but instead with a whole 
class of problems. Correspondingly, we propose that the 
overall problem be treated as a sequence of test prob- 
lems with increasingly more sophisticated treatment of 
the cavity boundary. Depending on the physical situa- 
tion, this sequence may also represent various stages in 
the time evolution of the system. 

For example, one can first consider the case where the 
cavity walls are rigid and have a fixed (e.g., spherical) 
shape (see § 12. ip . This may represent the early stages 



of the system's evolution. Next, one can assume that 
the shape and the size of the cavity are not fixed but in- 
stead are governed by the pressure balance between the 
electromagnetic stress inside the cavity and a constant 
external pressure outside (see § 12. 3p . This makes the 
set-up similar to Lynden-Bell's (1996) magnetic tower 
model. Thirdly, one can consider a pulsar in a fully dy- 
namic environment of a collapsing star. The set-up of the 
latter problem is essentially similar to that considered by 
Ostriker & Gunn (1971). 

All three versions of our pulsar-in-a-cavity problem are 
basic physics problems that ought to be solved if we are 
ever to understand how a millisecond magnetar works in- 
side a collapsing star. Our understanding of these prob- 
lems will benefit from rigorous mathematical analysis, 
but ultimately will most likely be achieved with the help 
of numerical simulations that are now becoming feasible. 
Whereas the first two problems represent perfect targets 
for relativistic force-free simulations, the third problem 
will most likely require a full relativistic MHD simula- 
tion. 

To set the stage for future numerical studies, and to be 
able to interpret their results, it is useful to get some ba- 
sic qualitative understanding of the problem. Therefore, 
in this section we will sketch what we think is a plausible 
physical picture of the system's evolution and how it re- 
lates to our magnetic tower model for GRBs (Uzdensky 
& MacFadyen 2006). 

In order to gain a more complete understanding of the 
interaction between the central magnetar and the sur- 
rounding stellar material, a full magnetohydrodynamic 
(MHD) description that includes plasma pressure and 
inertial effects will eventually be required. Of particu- 
lar interest would be the confinement of the expanding 
magnetosphere by the surrounding plasma and the dy- 
namical response of the star to the expanding magneto- 
sphere at its center. The full-MHD approach is especially 
relevant if there is a strong wind driven off the PNS by 
neutrinos and/or by the magneto-centrifugal mechanism, 
as considered by Thompson et al. (2004) and by Buc- 
ciantini et al. (2006). A useful simplification may come 
from noting that the main difference between the dense- 
plasma case and the relativistic force-free case is just 
the difference between the Alfven speed and the speed 
of light (J. Ostriker, private communication). Then, the 
MHD case may be treated similarly to the relativistic 
force-free case, but with the light cylinder replaced by a 
smaller Alfven surface. 

For simplicity, however, in this paper we shall restrict 
ourselves to the force-free case. That is, we shall assume 
that the plasma density inside the cavity is so low that 
electromagnetic forces dominate the dynamics almost ev- 
erywhere inside the cavity (but outside the neutron star 
of course) . The only exception is the part of the equato- 
rial plane outside the light cylinder, where plasma inertia 
needs to be taken into account (see below). While not 
realistic, given the large plasma densities present in the 
center of a massive star, the force-free description may 
nonetheless reflect some essential features of the full so- 
lution. It is of relevance especially for late phases of the 
evolution when the magnetic field outside the neutron 
star has been amplified to large values. 

As we have already mentioned, we shall also assume 
that the plasma inside the cavity can be accurately rep- 
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resented by an infinitely conducting fluid. We expect 
this key assumption to be well justified throughout most 
of the cavity, owing to the very large plasma densities 
and temperatures. Indeed, the high plasma density en- 
sures that the plasma (including photons) is highly col- 
lisional and hence is well described by resistive MHD; 
this means that the resistivity due to particle-particle or 
photon-particle collisions dominates over all other non- 
ideal terms in generalized Ohm's law. On the other hand, 
because of the very high plasma temperature, the resis- 
tivity is actually quite small, i.e., the magnetic Reynolds 
number is very high. All this makes ideal MHD a good 
approximation in the environment of a collapsing star 
(see Uzdensky & MacFadyen 2006 for more discussion). 
At the same time, we do acknowledge that this assump- 
tion may break down in some special regions, in partic- 
ular, inside the equatorial plasma strip (see below) and 
at the cavity boundary, where various fluid instabilities 
may lead to enhanced turbulent energy dissipation. In 
any case, however, we expect ideal MHD to be much 
better justified inside a collapsing star than the force- 
free assumption. For this reason, in this paper we shall 
ignore any finite-resistivity effects, leaving them for a 
future study. In addition, because of the very high den- 
sity, the plasma inside the cavity is completely optically 
thick to electromagnetic radiation and so the photons 
are tightly coupled to the gas. What this means is that 
there is no radiative (apart from possible neutrino cool- 
ing which we ignore in this basic problem) cooling in our 
system. Therefore, all the energy that is extracted from 
the pulsar stays inside the cavity, as long as we don't 
allow the cavity to expand. 

Throughout most of this discussion we shall ignore all 
numerical factors, e.g., 4ir, etc. Also, we shall assume 
that the magnetar rotation rate f2* stays approximately 
constant on the timescales under consideration. How- 
ever, ultimately one will have to consider the effect of 
decreasing rotation rate as the pulsar slows down. 

Finally, we would like to stress that there are impor- 
tant differences between the problem of an isolated pul- 
sar magnetosphere and our pulsar-in-a-cavity problem. 
In particular, in the isolated pulsar case one usually 
seeks a steady state (although perhaps employing time- 
dependent simulations to achieve it). In the pulsar-in-a- 
cavity case, on the other hand, we don't expect a station- 
ary solution; the problem is intrinsically time-dependent 
and it is the time evolution of the system that is of 
particular interest. In addition, it is believed that the 
wind of a normal isolated pulsar crosses the fast magne- 
tosonic surface somewhere far beyond the light cylinder 
and then reaches the termination shock. This is impor- 
tant because it implies that the inner pulsar magneto- 
sphere is causally disconnected from the outside; in par- 
ticular, the inner magnetosphere's structure and the pul- 
sar spin-down power cannot be influenced by the bound- 
ary conditions at very large distances. In sharp contrast, 
our case of a magnetosphere enclosed within a finite-size 
cavity is qualitatively different, because we now lack that 
huge separation of radial scales. In practical terms, this 
means that we draw a dividing line between the isolated 
pulsar magnetosphere and the bounded pulsar magne- 
tosphere based on the presence or absence of the fast 
magnetosonc surface inside the cavity. In particular, in 
our present study we are interested in the case of a cavity 
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formed inside the stalled supernova shock (see § [3]) . Its 
radius may be about 100-200 km, i.e., only moderately 
larger than the light cylinder radius of a millisecond mag- 
netar (about 30 km). Then, the plasma outflow just may 
not have enough range to reach the fast magnetosonic 
surface. As a result, our bounded pulsar magnetosphere 
always remains in causal contact with the outer bound- 
ary. Correspondingly, the inner structure of the magne- 
tosphere and the pulsar spin-down power is affected by 
the confining cavity. 

2.1. Pulsar in a Fixed Spherical Cavity 

We start with our problem I, in which the walls of 
the cavity are fixed. For definiteness, we take the cavity 
to be spherical in shape. The main results obtained in 
this section should also be approximately valid for the 
case of expanding (or contracting) cavity, as long as the 
expansion (contraction) speed is slow compared with the 
speed of light. 

Let us try to think physically about how the mag- 
netic field will evolve after the pulsar is spun-up instanta- 
neously at t = 0. In the Goldreich-Julian (1969) model 
for an isolated pulsar, the field lines extending beyond 
the light cylinder bend backwards and tend to become 
open. [Strictly speaking, the field lines actually always 
close, but very far away, in the so-called "boundary zone" 
(Goldreich & Julian 1969).] As long as the cavity bound- 
ary (the outer edge of the magnetosphere) lies outside 
the fast magnetosonic surface of the outflow, there is no 
feedback of this boundary on the inner magnetosphere. 
Then, the pulsar continuously spins down, losing its rota- 
tional energy and angular momentum to magnetic brak- 
ing by these effectively-open field lines. In our case, on 
the other hand, such an immediate field-line opening is 
not possible since the entire magnetosphere is contained 
inside the cavity of a finite size. This is one of the most 
critical differences between the isolated pulsar case and 
our case. 

Development of Differential Rotation 

One important point that one needs to take into ac- 
count is the establishment of differential rotation in the 
magnetosphere. This is nontrivial, since, by assumption, 
the magnetar rotates uniformly. However, as we will now 
show, the field lines that extend beyond the pulsar light 
cylinder nevertheless necessarily undergo differential ro- 
tation. As a result, these field lines are continuously 
twisted and hence toroidal flux is continuously injected 
into the cavity. 

To see how this comes about, let us consider a field 
line \E' (Fig. [2]) and compare the angular velocities at 
two points on this line: point A, where the field line at- 
taches to the pulsar, and point B, where it intersects 
the equator. Since this field line extends beyond the 
light cylinder, it cannot remain purely poloidal: toroidal 
field has to develop so that the plasma particles could 
slide backwards and out along the line, like beads on 
a wire. This toroidal field leads to a continuous brak- 
ing of the star so that there is an outward flux of angu- 
lar momentum and a Poynting flux of energy along the 
line. However, the toroidal magnetic field at point B 
has to be exactly zero because of the assumed reflec- 
tion symmetry with respect to the midplane. Therefore, 
the plasma can no longer slide toroidally; this means 
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Fig. 2. — Aligned pulsar inside an infinitely-conducting spher- 
ical cavity of radius Rq. After a time of order the light-crossing 
time Ro / c, the poloidal field lines outside the light cylinder expand 
somewhat but still remain confined within the cavity. Because the 
toroidal magnetic field has to vanish at the equatorial midplane 
due to reflection symmetry, the field lines there cannot corotatc 
with the star, Qg < c/Rg < f!„. As a result, differential rota- 
tion is established in both hemispheres, Af2 = Q* — Qg ~ f2» 
(for Rg ^> -RlcOj which leads to continuous generation of toroidal 
magnetic flux. 



that the toroidal velocity of the field line is equal to 
that of the plasma at this point. The angular momen- 
tum and rotational energy of the pulsar extracted by the 
magnetic field are partly accumulated and stored in the 
magnetic form and partly transferred to the equatorial 
plasma. Thus, the material at point B is continuously 
torqued by the magnetic field. Then, since the confin- 
ing wall prevents the material from moving out freely in 
the radial direction, the toroidal velocity of the plasma 
becomes closer and closer to the speed of light. How- 
ever, it can never exceed the speed of light; therefore, the 
plasma, and hence field-line, angular velocity at point B 
is bounded: fls — c/Rb = £!*i?Lc/-Rs- On the other 
hand, the angular velocity at point A is of course just 
the rotation rate of the pulsar: Qa = This means 
that the field line experiences differential rotation at a 
rate Ail = Q A - Sl B > - Rlc/Rb)- For field lines 
that cross the equator well outside the light cylinder, 
Rb Rlc, we then have Ail w f2». Thus, differen- 
tial rotation is established over a time-scale of order the 
light-crossing time across the cavity, to = Ro/c> fi" 1 . 

This differential rotation is important because it leads 
to a continuous toroidal stretching of the field lines and 
thus to a continuous injection of toroidal magnetic flux 
(of opposite signs) into the upper and lower hemispheres. 
Since all this toroidal flux has to be contained within a 
cavity of fixed size, the toroidal magnetic field at any 
given point grows, roughly speaking, linearly in time. 
This is in contrast with the pulsar in a free space, where, 
within a sphere of any given radius, a steady state is 
established on the time scale of order the light travel 
time across this radius. One worry that one might have 
in our case is the possibility of the kink instability as the 



magnetic field in the cavity becomes highly wound up. 
We address this issue in more detail in § 14.11 

Finally, the toroidal field reverses sharply across the 
equator, so there is a non-force-free equatorial current 
sheet that carries the radial return current back to (or 
from) the neutron star (see below). 



Now let us try to estimate the toroidal field evolution 
and distribution inside the cavity on long time scales 
(t = Nto, where N 3> 1, and to = Rq/c is the light 
crossing time across the cavity) and at distances much 
larger than the light cylinder radius. 

As we discussed above, because of the differential ro- 
tation, the bounded magnetosphere cannot be station- 
ary: toroidal magnetic flux is constantly being injected 
into a finite volume. Hence, the toroidal field strength 
continuously increases, whereas the poloidal magnetic 
field does not. The poloidal electric field, E po \, may 
become much larger than J5 po i but in any case cannot 
exceed the value B po iil^Ro/c — B po \Ro/ Rlc- Thus, af- 
ter several light-crossing times the magnetosphere out- 
side the light cylinder becomes toroidal-field dominated: 
B<p ^> E po i, B po i. 

Next, even though the configuration is time-dependent, 
after many light-crossing times the evolution slows down. 
Indeed, the poloidal field structure readjusts (e.g., in re- 
sponse to a change in the toroidal field strength) on a 
time scale of order the fast-magnetosonic crossing time 
across the cavity; for a force-free plasma this coincides 
with to = Rq/c. Since the toroidal flux grows linearly 
in time, the relative change in the toroidal field strength 
over At ~ to becomes small (of order TV" 1 ) at late times, 
t = Nto, N 3> 1. An approximate force- free equilibrium 
is then established separately in each of the two hemi- 
spheres, with poloidal current being approximately con- 
stant on poloidal flux surfaces: I ~ The magnetic 
field structure in such an equilibrium is governed by the 
relativistic force-free Grad-Shafranov equation (aka the 
pulsar equation). In the limit where the toroidal mag- 
netic field totally dominates the dynamics, this equation 
reduces to II' (Si?) = 0, so that the poloidal current func- 
tion becomes independent of , 5: 1(^6) = Io = const. This 
corresponds to the vacuum field produced by a singu- 
lar line current Io (which grows linearly in time) flowing 
along the rotation axis. The toroidal magnetic field is 
B$(t,R,Z) = I Q (t)/R, i.e., B$ — constant on cylinders, 
and the equilibrium can be described as the balance be- 
tween the toroidal field tension and pressure. In other 
words, the j x B force becomes relatively small inside 
the cavity, because the poloidal current becomes spatially 
separated from the toroidal magnetic field: it flows out of 
the pulsar along the axis (in both hemispheres) , then as a 
surface current along the cavity walls, and finally returns 
to the pulsar along the non-force-free equatorial current 
sheet. The bulk of the magnetosphere is thus almost 
current-free. In this regard, the electric-current struc- 
ture of the cavity is similar to that of the magnetic bub- 
ble considered by Lyutikov & Blandford (2002, 2003) in 
their model for Poynting-flux dominated GRB outflows 
(although we apply our model deep inside the collaps- 
ing star, that is, on different spatial and temporal scales 
compared with their model). 

Let us now estimate the magnitude of the poloidal line 
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current Io(t) and hence the characteristic strength of the 
toroidal field in the cavity. We shall express magnetic 
quantities characterizing the field in the cavity in terms 
of the total poloidal magnetic flux that extends beyond 
the light cylinder, which we shall call ^o- Up to a factor 
of order unity, this flux can be estimated from the pure 
dipole magnetic field, i.e., 



*0 ~ *dipole(-RLc) — 



Rlc 



(1) 



This estimate is justified because inside the light cylin- 
der the poloidal field remains close to dipole. Moreover, 
even in the extreme case of an unbounded, isolated pul- 
sar magnetosphcrc, in which the field is completely open 
outside the light cylinder, the poloidal flux crossing the 
light cylinder differs from the dipole formula only by a 
small amount (e.g., Contopoulos, Kazanas & Fendt 1999; 
Komissarov 2006; McKinney 2006b; Spitkovsky 2006). 
Thus, this estimate should be quite good in our case as 
well. 

Now, what is the characteristic poloidal magnetic field 
strength in the cavity at distances r ~ Rq ? Here the 



dipole formula (B po i 



), describing a fully-closed 



non-rotating field, and the split-monopole formula, de- 
scribing the fully-open magnetosphere of an isolated pul- 
sar, differ. In our case, all the field lines are closed, i.e., 
they intersect the equator within Ro, so one might think 
that the dipole-field estimate should be more applicable. 
However, as we show below, almost all of the field lines 
crossing the light cylinder actually intersect the equator 
in a narrow strip near the outer wall; therefore, the char- 
acteristic poloidal field at distances of order i? from the 
center and off the equatorial plane should be estimated 
as 

Bpoi ~ B = —J ■ (2) 
Ho 

For given by equation (fl]) , this estimate gives a value 
-Bpoi ~ B* (RI/RqRlc), which is by a factor R /Rlc 
larger than a pure dipole field at these distances. 

Now let us estimate the poloidal current and the 
toroidal magnetic field. In general, the poloidal current 
flowing through a region enclosed by an axisymmetric 
flux surface '3/ can be calculated by following the shape 
of a field line corresponding to "J: 



J(*,t) = Ant 



dl 



pol 



B po ii? 2 (; po i) 



(3) 



where Z po i is the path-length along the poloidal field. The 
main contribution to the integral comes from large dis- 
tances, R ~ Rq and thus the integral can be estimated as 
being of order i? / v E'o- Then, since Af2 ~ f2* = c/R^c, 
the axial poloidal current can be estimated as 



I (t) ~ Q*t 



*n t 



Rq i?LC to 



(4) 



Thus we see that for t 3> to the poloidal current becomes 
much stronger than the typical poloidal current in the 
unbounded pulsar magnetosphere (/ ~ ^o/Rlc)- Using 
the estimate ([1]) for V^o, we can express Iq as 



Io{t) ~ 5* 



Rt t 



(5) 



Correspondingly, the characteristic toroidal magnetic 
field at distances of order Rq is 



B<t>{Ro) 



Ro 



BoQ*t , 



(6) 



which is similar to the estimate presented by Kardashev 
(1970) for the toroidal field of a pulsar inside an expand- 
ing supernova cavity. We see that, after many light- 
crossing times across the cavity, B^Ro) becomes much 
larger than the toroidal field of an isolated pulsar at these 
distances [B™ utcd ~ f /(i? i?Lc) = B q (Ro/Rlc) = 
B n.*t < B tt*i\. 

Finally, we would like to remark on how to deter- 
mine the structure of the poloidal field, \&(r, 8). Usu- 
ally, when studying steady-state axisymmetric magneto- 
spheres, one uses an iterative procedure (e.g., Contopou- 
los et al. 1999). First, one makes a guess for the poloidal 
current J(\&), plugs it into the Grad-Shafranov equation, 
and solves this equation for ^(r, 8). Then one uses equa- 
tion ((3|) to determine the new function J(\P) and repeats 
the steps until the procedure converges. In our case, 
however, this approach does not appear to be feasible, 
since to lowest order the Grad-Shafranov equation sim- 
ply gives I(^) = Jo = const. We therefore advocate for 
an inverted approach where the poloidal flux function is 
determined from equation ([3]) . How to realize such an ap- 
proach in practice is not clear. One thing to note though 
is that this calculation should depend on ty(R,Z = 0) 
as a boundary condition, and this has to be determined 
from considering the redistribution of the poloidal flux 
across the equatorial midplane. This issue is discussed 
in the next subsection. 

Centrifugal Force in the Equatorial Plane 

As we noted above, the magnetosphere outside the 
pulsar light cylinder cannot be entirely force-free. Be- 
cause the toroidal magnetic field reverses across the equa- 
tor (due to the assumed reflection symmetry), the mag- 
netic field tension continuously accelerates the equatorial 
plasma in the toroidal direction. Correspondingly, this 
tension force performs mechanical work on the equato- 
rial plasma and so a certain part of the rotational energy 
extracted from the pulsar by the magnetic field is de- 
posited in the equatorial plane (the rest is stored in the 
bulk of the cavity as the toroidal magnetic field energy). 
Since the plasma in the equatorial plane rotates ultra- 
relativistically, the added energy leads to an increase in 
the relativistic "mass" of the plasma, Am ~ t 2 . An im- 
portant consequence is that this relativistically-rotating 
massive equatorial sheet experiences an outward radial 
centrifugal force, -Fccnt- This force cannot be balanced by 
the toroidal magnetic field because the latter is zero at 
the equator. Consequently, the equatorial plasma moves 
towards the wall and compresses poloidal magnetic field, 
until finally the centrifugal force is balanced by the j x B 
force due to the non-force-free part of the toroidal cur- 
rent j,/,. Thus, the poloidal magnetic flux in the equa- 
torial plane outside the light cylinder is pushed against 
the wall and is strongly concentrated in a narrow band 

4 The contribution from the the force-free part, j^[z = 0) = 
p e v^(z = 0) is exactly canceled by the radial electric force, p e E r , 
provided that the ideal-MHD condition cE = v X B holds in the 
equatorial strip. 
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Fig. 3. — At late times, the poloidal magnetic field is pressed 
against the wall by the centrifugal force of the rotating massive 
equatorial sheet. 



of ever-decreasing width d(t) <C Ro near the wall (see 
Fig. [3]). Because of this effect, we can expect nearly all 
the poloidal flux ^f that extends beyond the light cylin- 
der to cross the equator inside this strip, i.e., at cylin- 
drical radii R ~ Rq. At the same time, in the magneto- 
sphere above and below the equatorial plane, the poloidal 
field lines that emanate from this band have to fan out 
because they have to fill the cavity volume. Thus, the 
characteristic poloidal magnetic field in the cavity is of 
the order B = \& / Rq (see eqn. [5J and is much weaker 
(by a factor of d/R ) than in the equatorial strip. 

Let us assess the centrifugal force quantitatively. The 
total torque exerted on the massive equatorial strip by 
the magnetic field is given by r(t) = J I(^,t)d^ ~ 
Io(t)^o- Since the toroidal velocity of this strip is close 
to the speed of light, the total work per unit time due 
to this torque (i.e., the total Poynting flux that arrives 
at the strip) is P s trip — tc/Rq — B^c^/q. This power 
goes into accelerating the rotation of plasma in the strip, 
and some part of it may in principle dissipated into heat. 
Since the rotation here is already ultra-relativistic, the 
result of this acceleration is an increase in the rotation 
and/or thermal 7-factors, i.e., of the relativistic mass m 
of the plasma in the strip: d/dt(mc 2 ) = P s tri P - As a 
result, the relativistic mass grows with time as 



m(t)c 2 



Ro Ro 



Rlc 



Ro 



Rlc Ro Ro 



^ ~ Bl{t)Rl , 



that is, the plasma energy in the equatorial strip always 
remains small compared with the energy B^(t)R^ stored 
in the toroidal magnetic field at these distances. The 
centrifugal force acting on the equatorial strip can be 
estimated as 



t(*) = 



m(t)c 2 
Ro 



(8) 

We see that this force grows quadratically with time, just 
as the toroidal field pressure, but always remains small 
(by a factor of Rlc/Rq <C 1) compared with the overall 



horizontal force exerted on the side wall by the toroidal 
field. 

A detailed analysis of the internal structure of the mas- 
sive equatorial plasma strip, including its vertical struc- 
ture, lies beyond the scope of this paper. However, we 
present here a simple estimate for the Lorentz factor due 
to rotation, 7 rot , in terms of the strip width d and half- 
thickness h. This estimate is derived under a certain very 
restrictive set of assumptions and serves for illustration 
only. 

Let us consider the vertical force balance inside the 
strip in the co-rotating frame, and let us neglect the con- 
tribution from electric force for simplicity. Then the 
toroidal magnetic field pressure outside the strip has 
to be balanced by the plasma pressure inside: p co — 
B 2 /8ir. Next, let us make the assumption that the 
plasma in the strip is a light relativistic fluid with the 
adiabatic index 4/3. Then, the co-moving energy den- 
sity is pcoC 2 = 3p co . On the other hand, the total plasma 
energy mc 2 inside an annular strip of radius Ro, width d, 
and thickness 2h can be written in the lab frame as 
mc 2 — AirRftdh p co c 2 j 2 ot . By combining all these expres- 
sions with the equation ([7]) for mc 2 , we find 



7rot 



2 R Rlc 

3 dh 



(9) 



On the other hand, it may be possible that a significant 
amount of plasma accumulates in the equatorial strip 
or that the plasma there is strongly compressed by the 
toroidal field pressure. Then, the baryon number den- 
sity rib may become so large that the co-moving energy 
density is dominated by the non-relativistic component, 



i.e., by the baryon rest-mass, p c 



n btCO m p c* > p 

C 



B^/8tt. Since n 6 , co = 7 rot 
true can be written as 



I nt,, the condition that this is 



7rot > ° = 



B 2 



^n bmp ~ 2 - ^ 

Provided that we are in this regime, the total plasma 
energy in the strip is dominated by the kinetic energy 
of the baryons: mc 2 = 4irRodhj m tnbm p c 2 . Then, using 
equation ([7]), we get 



7rot 



RoRlc 
dh 



^strip • 



(ii) 



By substituting this expression into the condition (|10p . 
we see that the co-moving energy density is dominated 
by the rest-mass of the baryons only when 



dh 



RlcRo 



< 1 



Correspondingly, we have 



7rot < 



RoRlc 



dh 



(12) 



(13) 



Magnetic Spin-down Power of a Pulsar in a Fixed 
Cavity 

Another extremely important point is that the rate at 
which the magnetic field in a bounded magnetosphere ex- 
tracts rotational energy from the central rotating conduc- 
tor actually grows with time. This is because the mag- 
netic torque per unit area is proportional to the toroidal 
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field at the conductor's surface and the latter grows lin- 
early with time. Thus, the magnetic power generated by 
a spinning pulsar inside a cavity increases linearly with 
time as long as the cavity does not expand (or expands 
slowly) and the spin rate of the pulsar stays constant. 
We can estimate the spin-down power as 

P(t) = J(t)* ^* = flj* IT ~ Elated £ , (14) 

Ho Ho 

where Pi SO iatcd ~ B^R^Q^/c 3 is the spin-down power of 
an isolated, unbounded pulsar. As we see, after many 
light-crossing times, the power of a pulsar-in-a-cavity 
greatly exceeds that of a classical isolated pulsar. This is 
our answer to the apparent paradox raised by Lyutikov 
(2006). 

We thus emphasize that the energy extraction from a 
magnetar-in-a-cavity can be a run-away process. This 
is because the twisting of a magnetic field confined by 
an external boundary results in an increase in the field's 
strength at the light cylinder and hence in a growing rate 
of energy extraction from the magnetar. 

This effect can be attributed to a positive feedback 
that exists between the energy that has been already ex- 
tracted from the pulsar, and the strength of the agent 
that extracts the energy (the toroidal magnetic field). 
Namely, most of the extracted energy is stored in the 
toroidal magnetic field, and since the volume occupied 
by this field is kept finite, the toroidal field strength in- 
creases with time. Because the magnetosphere remains 
in a quasi-equilibrium, the toroidal field constantly read- 
justs everywhere, including within the light cylinder. In 
other words, because the system is not hyperbolic but 
elliptic, the inner magnetosphere feels the presence of 
the outer confining wall. In particular, the toroidal field 
at the very surface of the pulsar increases linearly with 
time, and hence so does the spin-down torque exerted by 
the magnetic field on the pulsar. This picture is similar 
to what is happening in the combustion chamber of a 
rocket, for example. In that case, the gas temperature 
and pressure increase as the chemical energy of the fuel 
is released in the combustion process. At the same time, 
the rate at which fuel burning occurs increases with an 
increase in the ambient temperature. As a result, rapid 
and efficient burning demands high pressure and hence 
a strong confining chamber capable of withstanding this 
pressure. Similarly, in our case of a pulsar placed inside 
a cavity, the presence of strong cavity walls leads to an 
increased energy extraction rate from the pulsar. 

In a realistic situation, this steady power growth might 
not last indefinitely. It may be limited, for example, 
by the development of the kink instability, which would 
result in the conversion of the toroidal flux to poloidal 
flux and to partial dissipation of magnetic energy (see 
§ 14.11 for more discussion) . 

2.2. Hoop-stress collimation: contrast with the isolated 

pulsar 

The toroidal field generated by the differential rotation 
exerts a constantly-growing pressure on the cavity walls. 
If we now relax the assumption that the walls are fixed 
and allow them to move, this pressure will make the cav- 
ity inflate. We then want to understand how rapidly such 
inflation will proceed and whether it will be isotropic or, 



say, collimated along the axis. We discuss the collimation 
issue in this subsection. 

Generally speaking, since the toroidal field pressure in 
the lateral direction is partly negated by the field's ten- 
sion (the hoop stress), which has no vertical component, 
one may expect the resulting expansion to be predomi- 
nantly vertical. However, notice that here we are inter- 
ested in a situation where the (differential) rotation is 
relativistic: AflRo ~ fl*Ro S> c. On the other hand, 
Lynden-Bell's (1996) magnetic tower model, for exam- 
ple, was developed for the non-relativistic regime. It is 
well-known that hoop-stress collimation is not a trivial 
issue in the relativistic case. Thus, it is not immediately 
obvious that the hoop-stress collimation mechanism can 
be applied to the pulsar-in-a-cavity scenario considered 
in this paper. We therefore would like to discuss this 
issue in some detail here. 

At first, one might think that there should be no prob- 
lem collimating the outflow: the magnetic field is pre- 
dominantly toroidal even without differential rotation. 
And it is the toroidal field's hoop stress that is usu- 
ally credited for collimating astrophysical jets. However, 
as is well known, hoop-stress collimation does not work 
as well when applied to ultra-relativistic magnetically- 
dominated outflows, as it does in the non-relativistic 
case. The quintessential example of this lack of colli- 
mation is the isolated-pulsar wind inside the termina- 
tion shock. The basic reason for this difficulty is the 
decollimating force due to the poloidal electric field, 
E po \. Indeed, in the case of an unbounded relativistic 
uniformly-rotating force- free magnetosphere (e.g., an iso- 
lated aligned pulsar magnetosphere) in a steady state, 
the poloidal electric and toroidal magnetic fields have to 
be nearly equal in strength at large distances from the 
central axis (Goldreich & Julian 1969). Importantly, it 
turns out that this balance can be realized in an uncol- 
limated, quasi-spherical poloidal magnetic field configu- 
ration; an excellent example of this is Michel's (1973) 
split-monopole solution. A rough argument explaining 
this lack of hoop-stress collimation in the relativistic- 
rotation case goes as follows. Let us consider an uncolli- 
mated field configuration; the poloidal magnetic field is 
open outside the light cylinder and has a split-monopole 
geometry, i.e., drops off with distance as r~ 2 . In a steady 
state, the poloidal electric field is E po \ — B po \ R/ Rlc 
where R is the cylindrical radius. It therefore drops off 
along radial rays as r _1 . But the toroidal magnetic field 
also drops off as r _1 . Moreover, at the light cylinder E po \ 
and are comparable: E po \ — B po \ ~ B^. Since out- 
side the light cylinder they both decrease as the same 
power of r, they remain comparable to each other (both 
being much larger than B pQ \) at large distances. More- 
over, as Goldreich & Julian (1969) showed, E po \ and B$ 
even become equal asymptotically as r — > oo. The bot- 
tom line is that a quasi-spherical relativistic force-free 
equilibrium can be established as a balance between the 
collimating pinch force (the sum of the toroidal magnetic 
field pressure and its tension) and the opposing electric 
force. Hoop-stress collimation is suppressed as a result 
of this balance. 

Now, in the case of a rotating magnetosphere enclosed 
inside a rigid cavity of a fixed radius Rq > i?LCj the 
situation is different and hoop-stress collimation can in 
fact work. Indeed, as we showed above, after many 
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light-crossing times it ^> Ro/c), the toroidal magnetic 
field filling the cavity becomes stronger than both B po i 
and -Epoi, in contrast to the isolated pulsar case. More- 
over, this toroidal field is distributed nonuniformly; it 
is basically inversely proportional to the cylindrical ra- 
dius. Correspondingly, the stress exerted by the toroidal 
magnetic field on the cavity walls is also nonuniform: 
the magnetic pressure pushing vertically against the top 
and bottom walls is much higher that the lateral mag- 
netic pressure acting on the side walls. Therefore, if 
we now allow the cavity to expand under this pressure, 
we expect any subsequent expansion to be mostly ver- 
tical (see Fig. SJ, at least as long as the expansion ve- 
locity is slow compared with the speed of light. Then 
we effectively find ourselves in a situation similar to the 
non-relativistic magnetic tower proposed by Lynden-Bell 
(1996). We therefore envision that the eventual, long- 
term result of this process will be the creation of a pair of 
oppositely-directed magnetic towers (Uzdensky & Mac- 
Fadyen 2006). The interaction of the expanding towers 
with the surrounding stellar envelope aids in their con- 
finement, similarly to jet collimation seen in hydrody- 
namical simulations of the collapsar model (MacFadyen 
& Woosley 1999; Aloy et al. 2001; MacFadyen, Woosley 
& Heger 2001; Zhang, Woosley & MacFadyen 2003). In 
the scenario considered in the present paper, these tow- 
ers are driven not by a differentially-rotating disk, but by 
a rapidly-rotating magnetar. This suggests that consid- 
ering the magnetosphere of a pulsar inside a cylindrical, 
as opposed to spherical, cavity may represent yet an- 
other interesting and important problem for future study 
(see § S3). 

An important element in the above discussion is the 
fact that the electric field is small compared with the 
toroidal magnetic field. This is directly related to the 
fact that the toroidal field is generated not as a part of 
an outgoing large-scale electromagnetic wave driven by 
the pulsar rotation, but as a result of differential rotation. 
This observation points to the important role played by 
differential rotation (as opposed to uniform relativistic 
rotation) in collimating relativistic force-free outflows. 

2.3. Pulsar magnetosphere confined by a constant 
external pressure 

Let us now consider the case when the pulsar magne- 
tosphere is confined by some fixed and uniform external 
gas pressure, P G xt, instead of a cavity of fixed radius Rq. 
We are interested in this particular set-up because it is 
closest to that considered by Lynden-Bell in his original 
magnetic tower paper (Lynden-Bell 1996), and we here 
want to compare his non-relativistic disk model with a 
pulsar in a similar setting. 

Like Lynden-Bell, let us assume that the external pres- 
sure is weak compared with the magnetic field pressure 
B% I 8tt in the immediate vicinity of the rotating conduc- 
tor. Moreover, because we are interested in exploring rel- 
ativistic effects, we want our pulsar magnetosphere to be 
able to expand well beyond the light cylinder. Therefore, 
we shall also assume that the external pressure is small 
compared with the magnetic pressure of a pure dipolc 
field at the light cylinder: 87rP ex t < Bl ipole (R LC ) ~ 

BUR*/RlcY 
Let us imagine, as is frequently done in time-dependent 



t» 




FlG. 4. — Axisymmctric pulsar inside a cavity. After many light- 
crossing times, the magnetosphere becomes toroidal-field domi- 
nated. Because of hoop-stress, the magnetic stress on the cavity 
becomes strongly concentrated near the axis. This leads to a pre- 
dominantly vertical, collimated expansion; a magnetic tower forms. 



pulsar magnetosphere studies (e.g., Komissarov 2006; 
McKinney 2006b; Spitkovsky 2006), that we start with 
a non-rotating star with a dipole field and then spin it 
up suddenly at t = 0. The initial evolution of the mag- 
netic field is then similar to that of an isolated pulsar: 
the field lines that extend beyond the light cylinder start 
to wind up and expand radially at the speed of light, 
i.e., i?o — ct. This stage of uninhibited quasi-spherical 
expansion proceeds until the magnetic field pressure at 
the outer edge of the expanding magnetosphere becomes 
as small as the external gas pressure. In order to esti- 
mate when this happens, we need to evaluate the toroidal 
magnetic field pressure at R = Ro(t). The toroidal 
field changes with time because of two opposing fac- 
tors: continuing injection of the toroidal magnetic flux, 
x{t) ~ ^o^** = ^oRo{t)/RhC, and the the increasing 
volume of the cavity. The net result is that the toroidal 
field drops off according to 



B4R (t),t] 



x(t) 



R 2 (t) R hC Ro(t) Rhcct 



(15) 



Another way to obtain this estimate is to note that the 
main result of this free expansion is the establishment of 
the stationary isolated-pulsar magnetosphere inside the 
radius Ro(t). The toroidal magnetic field in an isolated 
pulsar magnetosphere scales as B^(r) ~ ^o/Rlct" (Gol- 
dreich & Julian 1969), which is equivalent to the above 
estimate. 

Eventually, the pressure of the toroidal magnetic field 
drops to a level where it becomes equal to the exter- 
nal gas pressure (there is also a comparable contribution 
from the electric field). This happens at time t = t cqi 
corresponding to the cavity radius reaching an equilib- 
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rium value 



R — rt = 

V87rP oxt i? LC 



(16) 



After this, the expansion continues, but changes its 
character: the lateral expansion slows down, and the ex- 
pansion becomes mostly vertical. Eventually, at t 3> t eq , 
a magnetic tower forms, similar to Lynden-Bell's (1996) 
tower. One important difference is that since the radius 
of the tower is much larger than i?LCj & proper analy- 
sis requires relativistic treatment, so Lynden-Bell's non- 
relativistic theory is not directly applicable. In partic- 
ular, we expect the vertical expansion of the tower to 
be relativistic. This can be seen from the following ar- 
gument. As the tower grows, its radius stays roughly 
constant, of order R eq , whereas its height increases lin- 
early with time, with the velocity Vt op - The continuously 
injected toroidal flux goes into filling the expanding vol- 
ume of the tower with toroidal magnetic field, so that, 
roughly speaking, 



B<l> Hop t R oq . 



(17) 



Assuming ~ \/87rP cx t, we therefore arrive at the es- 
timate 

Vtop~c. ( 18 ) 

This result can be understood naturally by noting that 
the problem has no mass or density parameter and so 
there is no characteristic velocity scale other than the 
speed of light c (scales like f2*i? e q are even larger than c). 

The toroidal magnetic field stays roughly constant dur- 
ing this stage, and so the poloidal current flowing through 
the tower is also constant and is of order 



_*0_ 



(19) 



the same as the poloidal current in the isolated pulsar 
case. Correspondingly, the magnetic luminosity (i.e., the 
spin-down power of the pulsar) stays at a constant level of 
order isolated- However, unlike the isolated pulsar case, 
this luminosity is not quasi-spherical, but is channeled 
predominantly in the vertical direction. 

Provided that the expansion of the tower is sub- 
magnetosonic, an approximate relativistic force-free 
equilibrium is established inside the tower (at least away 
from the top lid of the tower). As in the fixed-cavity 
case, the work done by the toroidal field's magnetic ten- 
sion on the equatorial current sheet goes into accelerating 
the equatorial plasma to ultra-relativistic velocities. The 
relativistic mass of this plasma, and hence also the ra- 
dial centrifugal force grow linearly with time, as does the 
overall magnetic pressure force on the outer wall (because 
of the steadily increasing height of the tower). 

3. MAGNETAR INSIDE A COLLAPSING STAR: AN 
OUTLINE OF THE GENERAL SCENARIO 

Previous studies of core-collapse supernovae (SNe) 
have shown that, when the core of a massive star 
collapses into a proto-neutron star, a bounce shock 
is launched back into the star (see the reviews by, 
e.g., Bethe & Wilson 1985; Woosley & Weaver 1986; 
Bethe 1990). However, as was also shown in these stud- 
ies, the shock quickly stalls at a radius of about 200 km. 
The explosion then enters a relatively long (~ 1 sec) 
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Fig. 5. — Stalled shock phase of core-collapse explosion. 



quasi-stationary phase (see Fig. [5]) . During this phase ac- 
creting material constantly moves through the shock and 
gets heated to very high temperatures. The shock looks 
stationary in the Eulerian frame and the shock jump con- 
dition can be viewed as a balance between the ram pres- 
sure of the infalling material, that tends to quench the 
shock, and the thermal pressure of the post-shocked gas, 
that is supported mostly by the continuous heating due 
to neutrino deposition in the dense plasma behind the 
shock. Gradually, both the neutrino luminosity and ac- 
cretion rate decline with time. Eventually, one of two 
things has to happen as an outcome of the competition 
between neutrinos and accretion. If neutrinos win, the 
shock engulfs the entire star and one gets a successful 
SN explosion. If they lose, the shock dies and the PNS 
gains mass beyond the critical mass and collapses into a 
black hole, which then subsequently swallows the rest of 
the star, without a SN. 
In our model, we add a third dynamical component 
- the magnetic field. The magnetic force is pushing 
out, helping the explosion, as is the thermal pressure 
of the neutrino-heated gas. These two outward forces 
are opposed by the accretion ram pressure. Our main 
idea is that, generally speaking, the two outward forces 
evolve differently with time, and thus the explosion may 
be a two-stage process. In particular, we suggest that 
the magnetic pressure force is unimportant during the 
stalled-shock phase that lasts a few hundreds of msec. 
However, we note that during this time the magnetar 
makes several hundred revolutions, resulting in a great 
amplification of the toroidal magnetic flux by the differ- 
ential rotation. [Of course, during this stage the field 
is not force-free, and the gas pressure and inertia are 
important.] Over time, however, both the neutrino en- 
ergy deposition and the accretion rate decline, whereas 
the toroidal magnetic field grows (see Fig. [6]). For ex- 
ample, assuming Rq = 3i?LC = 10-R* = 100 km, and 
-B* = 10 15 G, we see that the entire cavity is filled with 
3 • 10 14 G fields after about 100 turns (0.1 sec), cor- 
responding to the magnetic pressure of about 4 • 10 27 
erg/cm 3 . This is to be compared with the ram pressure 
of the infalling stellar material that tries to compress the 
magnetosphere. The simplest estimate of the ram pres- 
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Fig. 6. — Schematic time evolution of the main three forces 
responsible for the stallcd-shock force balance. 



sure at r — i?n is given by 

P ram ~ ¥™ c 8 • 10 27 M Ml' 2 R^crg cm" 3 , (20) 
47rrtQ 

where v s = (2GM/i? ) 1/2 ~ 5 • 10 9 M 1/2 R^' 2 cm/sec is 

the free-fall velocity at radius Rq, and M and M are the 
mass enclosed within radius Rq and the accretion rate 
at this radius, expressed in units of M Q and M Q /sec, 
respectively. Thus, after a few hundreds of millisec- 
onds, the magnetic pressure overtakes the rapidly de- 
creasing neutrino heating as the main driving force and 
re-energizes the stalled shock, leading to a successful ex- 
plosion. This scenario is consistent with the picture pre- 
sented by Akiyama et al. (2003) who demonstrate numer- 
ically the growth of the magnetic field on the 200 msec 
timescale up to about 10 15 G in the range of radii up to 
100 km. The overall outcome of scenario is also similar 
to that suggested by Bucciantini et al. (2006), although, 
because of the winding-up amplification, the magnetic 
field becomes dynamically important much sooner in our 
model. 

To summarize our picture, the ram pressure of the 
accreting material provides a nurturing womb in which 
the baby magnetic field grows, until it is finally strong 
enough to break out. Neutrino energy deposition plays 
an important role during this gestation period, as it 
provides the support that prevents the magnctosphcrc 
from being completely squashed and buried by the ac- 
creting gas. Finally, if the above picture is correct and 
the explosion does become magnetically-driven, then the 
hoop-stress mechanism makes it highly collimated, thus 
satisfying one of the key necessary conditions for GRB. 
Note that this jet is driven by the magnetar-level (i.e., 
~ 10 15 G) field and is therefore stronger and faster 
than the LeBlanc- Wilson (1970) jet that may have been 
launched a few seconds earlier, during the core-collapse 
process (Wheeler et al. 2000). 

4. DISCUSSION 
4.1. Effect of MHD Instabilities 



The physical picture presented in this paper, with its 
smooth coherent magnetic structure, is, of course, an ide- 
alization, necessary for obtaining a basic physical insight 
into the system's dynamics and for getting the main ideas 
across in the clearest possible way. The actual mag- 
netic field is likely to be different from such a simple 
system of nested axisymmetric flux surfaces. Instead, 
it may consist of many loops of different sizes and ori- 
entations. It may thus have a complex substructure on 
smaller scales, both temporal and spatial. This substruc- 
ture may arise naturally from the beginning, especially 
if the proto-magnetar's magnetic field is produced by a 
turbulent dynamo. On the other hand, it may also result 
from a nonlinear evolution of various MHD instabilities 
that may develop in the system. The effect of MHD in- 
stabilities in our highly twisted magnetosphere is one of 
greatest uncertainties in our model. This section is de- 
voted to the discussion of two such instabilities: kink and 
Rayleigh-Taylor . 

(i) As the confined magnetosphere is twisted up, it may 
become prone to a non-axisymmetric kink-\ike instabil- 
ity. This may happen both during the pulsar-in-a-cavity 
phase and during a later magnetic tower phase. 

The kink is probably the most dangerous instability 
in our scenario. In its nonlinear stage, it may lead to 
a significant disruption. Such a disruption, however, is 
not necessarily a bad thing: it is likely to be only tem- 
porary and the tower may be able to reform after be- 
ing disrupted, as is seen in laboratory experiments by 
Lebedev et al. (2005). The resulting non-steady evolu- 
tion may then provide a plausible mechanism for rapid 
variability seen in gamma-ray bursts. In addition, as 
a result of such disruption, a fraction of the toroidal 
magnetic field energy may be dissipated into thermal en- 
ergy (Eichler 1993; Begelman 1998). As was shown by 
Drenkhahn & Spruit (2002; see also Giannios & Spruit 
2005, 2006, 2007), this may contribute to the acceleration 
of the Poynting-flux dominated outflow and to powering 
the prompt gamma-ray emission at later times. Also, 
such kink-driven magnetic dissipation in the magneto- 
sphere may be seen as a manifestation of "coronal ac- 
tivity" that may modify (harden) the emitted neutrino 
spectrum (Ramirez-Ruiz & Socrates 2005). 

As far as we know, the stability of the pulsar-in-a- 
cavity has not yet been studied. However, several non- 
relativistic 3D MHD simulations (Kato et al. 2004b; 
Nakamura et al. 2007; Ciardi et al. 2007) have recently 
addressed the kink instability of magnetic towers (al- 
though not in the GRB context). They seem to indi- 
cate that during the first few rotation periods, a tower 
is stabilized by the surrounding high-pressure gas, but 
at later times a large-scale external kink does develop. 
As a result, the tower's overall shape becomes helical. 
This, however, does not immediately lead to the total 
disruption of the tower; even though the configuration 
is nonaxisymmetric, its main morphological features re- 
main similar to those in the axisymmetric case (Naka- 
mura et al. 2007) . Similar conclusions have been reached 
by Nakamura & Meier (2004) in their 3D-MHD study 
of Poynting-flux-dominated jets propagating through a 
stratified external medium. These authors found that 
the jet stability strongly depends on the background den- 
sity and pressure profiles along the jet. In particular, a 
steep external pressure gradient forestalls the instabil- 
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ity onset. When the instability does eventually develop, 
the resulting helical structures saturate and do not de- 
velop into full MHD turbulence. An important theoret- 
ical evidence supporting external pressure stabilization 
follows from Konigl & Choudhuri's (1985) analysis of a 
force-free magnetized jet confined by an external pres- 
sure. They argued that a non-axisymmetric helical equi- 
librium state becomes energetically favorable (conserving 
the total magnetic hclicity in the jet) only when the pres- 
sure drops below a certain critical value. If this happens 
and the external kink mode does become unstable, then 
this non-axisymmetric equilibrium can be interpreted as 
the end point of the non-linear development of the insta- 
bility. 

In addition to the above non-relativistic studies, a 
few first steps have recently been taken towards under- 
standing the stability of relativistic jets, in particular, in 
the framework of relativistic force-free electrodynamics 
(Gruzinov 1999; Tomimatsu et al. 2001). However, to the 
best of our knowledge, to date there have been no formal 
stability studies of relativistic magnetic towers or of con- 
fined pulsar magnetospheres. Such studies, both analyti- 
cal and numerical, are clearly needed. They may involve 
a linear perturbation analysis or a non-axisymmetric rel- 
ativistic MHD or force-free simulation. They would have 
to take into account several stabilizing effects. First, as 
Tomimatsu et al. (2001) found in their linear stability 
analysis of a narrow rotating relativistic force-free jet, 
rapid field-line rotation inhibits kink instability. Sec- 
ond, regarding the stability of a rapidly growing magnetic 
tower, we expect that the tower expansion should quickly 
transition to the relativistic regime, eventually reaching 
a very large 7-factor. Once this happens, the relativis- 
tic time delay may effectively stabilize the outflow (see 
Giannios & Spruit 2006) . This is because MHD instabil- 
ities grow on the local Alfven-crossing time in the fluid 
frame and hence much slower in the laboratory frame. 
As a result, even if instabilities are excited, they do not 
have enough time to develop before the break-out of the 
flow from the star. 

We would like to point out that, in this problem, we are 
actually interested not so much in the instability onset 
or its early linear development, but rather in its long- 
term (many rotation periods) nonlinear evolution and 
its overall effect on the magnetosphere. Such a long- 
term behavior is very poorly understood and needs to 
be investigated in the future. Therefore, here we can 
only provide a hypothetical discussion. In its nonlin- 
ear stage, the kink instability may lead to conversion of 
some of the toroidal magnetic flux to poloidal flux (in 
our geometry). Some of this new poloidal flux may be- 
come detached from the star via reconnection (which, in 
reality, may be strongly inhibited deep inside the collaps- 
ing star, see Uzdensky & MacFadyen 2006). This could 
lead, in principle, to the break up of the single coherent 
magnetosphere into a number of smaller spheromak-likc 
plasmoids (Fig. [7|. That is, instead of further twist- 
ing up of the entire magnetosphere or, at a later stage, 
lengthening of the magnetic tower, one would effectively 
get continuous injection of new plasmoids into the sys- 
tem. Hoop stress still works inside each of them, and so 
the overall dynamical effect may be the same as that of 
a single magnetosphere, at least qualitatively. The re- 
sulting multi-component structure of the outflow may be 




Fig. 7. — Magnetic tower may have a substructure represented 
by a train of many spheromak-like plasmoids. This situation may 
arise as a result of spatial and temporal intermittency at the base 
of the outflow and/or due to instabilities and reconnection. 



responsible for the observed intermittency in GRBs. It is 
a very interesting scenario that should to be considered 
in future research. 

An important consideration that then needs to be 
taken into account is the conservation of magnetic helic- 
ity Differential rotation leads to a continuous injection 
of helicity into the system (of opposite signs in the two 
hemispheres). The kink instability may convert toroidal 
flux to poloidal but it will not destroy the magnetic he- 
licity accumulated in the cavity. Thus, whatever the re- 
sulting configuration might be, it will have to be consis- 
tent with a growing amount of helicity. One may in fact 
imagine a cyclic process involving twisting up the mag- 
netosphere for several rotation periods, followed by flux 
conversion due to the kink instability, followed by recon- 
nection and the production and detachment of a plas- 
moid carrying the magnetic helicity (and some of the 
magnetic energy) injected during the given cycle. One 
can hypothesize that if this cyclical process is robust, 
then over time the bulk of the cavity may become filled 
with spheromak-like plasmoids. The picture then would 
be somewhat similar to that in Figure [7l except it would 
not have to be axisymmetric. Each of these plasmoids 
would have some net magnetic helicity and magnetic en- 
ergy and would be in a magnetostatic equilibrium con- 
figuration, confined laterally by the overall pressure of 
the neighboring plasmoids. It may also contain thermal 
energy produced as a result of reconnection during plas- 
moid creation. Since the magnetic field is closed within 
each plasmoid, each plasmoid is not magnetically con- 
nected to the star and thus the magnetic field inside of 
it is not subject to any additional twisting. Helicity then 
stays constant within each plasmoid. As the number 
of such plasmoids grows with time, they occupy larger 
and larger fraction of the cavity volume. Correspond- 
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ingly, the part of the volume that is directly connected 
to the rotating neutron star shrinks with time. As long 
as some part of the neutron star field lines extends be- 
yond the light cylinder, differential rotation continues on 
these field lines, resulting in more twisting and more gen- 
eration of plasmoids. At some point, however, the part 
of the magnetosphere that is directly connected to the 
neutron star — the NS magnetosphere proper — will be 
squeezed by the surrounding plasmoids to such a degree 
that it will be confined entirely inside the light cylin- 
der. If this happens, twisting will stop, the neutron star's 
proper magnetosphere will be corotating with the neu- 
tron star and magnetic energy extraction will cease, at 
least if the boundary between the proper magnetosphere 
and the surrounding plasmoids is axisymmetric. If this 
boundary is not axisymmetric, then energy and angular 
momentum will continue to be extracted at some level 
through a process akin to the magnetic propeller effect 
(Illarionov & Sunyaev 1975). On the other hand, it is 
possible that magnetic reconnection will restore the link 
between the magnetosphere proper and the outer plas- 
moids, even if temporarily. If this happens, the situation 
will become more complicated; the spin-down torques 
will again be modified. In fact, as was pointed out to 
us by the referee of this paper, this resulting situation 
may become somewhat analogous to the case of a pulsar 
magnetosphere compressed by a strong relativistic wind 
of another pulsar, such as in the case of the double pul- 
sar PSR J0737-3039. Such magnetosphere-wind interac- 
tion, along with the resulting pulsar spin-down torque, 
was considered by Lyutikov (2004) and by Arons et al. 
(2005). As they point out, the torque due to the recon- 
nected field lines can become much larger than the usual 
spin-down torque of an isolated pulsar. 

It is also interesting to make the following comment. 
An unbounded relativistic force-free outflow driven by 
a rotating conductor is expected to be stable. On the 
other hand, a closed confined magnetosphere with field 
lines subject to differential rotation, such as our pulsar- 
in-a-cavity problem or a magnetic tower, may be kink- 
unstable. At the same time, as we discussed in Sec. 12.21 
the outflow is uncollimated in the first case but is colli- 
mated in the second case. This suggests that there may 
be a deep connection between stability and lack of colli- 
mation for axisymmetric relativistic force- free flows. 

To sum up, the kink, and especially its nonlinear out- 
come, is a serious issue that needs to be addresses in 
the future. Axisymmetric mode (i.e., sausage instabil- 
ity) also needs to be investigated. 

{it) Another important process that may potentially 
plague the development of coherent magnetic structures 
inside collapsing stars is the development of the Rayleigh- 
Taylor instability (or, rather, its magnetic counter-part, 
the Kruskal-Schwarzschild instability). This instability 
is expected to affect purely hydrodynamic fireball models 
as well; a strong magnetic field will suppress it somewhat, 
although probably not completely. As Wheeler ct al. 
(2000) and Arons (2003) have pointed out, the Rayleigh- 
Taylor instability at the interface between the lightweight 
relativistic fluid (electromagnetic field and hot relativis- 
tic plasma) and the overlying colder, denser stellar ma- 
terial may cause splitting of a coherent magnetic struc- 
ture into several separate strands interlaced with stellar 
matter. In the Arons (2003) model, the stellar enve- 



lope is quickly "shredded" by the nonlinear Rayleigh- 
Taylor "fingers". In effect, these fingers play a role of 
evacuated channels that allow the electromagnetic rela- 
tivistic plasma energy produced near the central engine 
to escape through the star. Arons further argued that 
these channels suffer only a relatively small amount of 
mixing with the non-relativistic stellar material due to 
the Kelvin-Hclmholtz instability. In light of this work, 
we cannot rule out the possibility that our magnetic 
cavity and/or the subsequently-formed magnetic towers 
may also suffer from fragmentation into several Rayleigh- 
Taylor "fingers". However, initially small-scale fingers 
quickly merge with each other to form a small number of 
large ones in the nonlinear stage. Therefore, we do not 
expect strong mixing of the baryons from the stellar enve- 
lope into the magnetosphere. The exact geometry of the 
outflow may change and a strong time-variability may 
develop, but, overall, we expect the outflow to survive. 
More research is needed in order to assess the implica- 
tions of this instability for our scenario. 

Finally, we would like to reiterate that a proper treat- 
ment of these problems requires a time-dependent 3D 
relativistic force-free or full (preferably relativistic) MHD 
analysis and simulations (see Sec. 14. 6|) . 

4.2. Nickel Production 

A central issue for the central engine of long-duration 
GRBs is the required production of 56 Ni. The supernovae 
that have been observed to accompany long duration 
GRBs (SN-GRBs) are classified as Type Ibc (SNe Ibc; 
see, e.g., Soderberg 2006; Kaneko et al. 2007). Modeling 
of the optical light curves of SNe Ibc requires the pres- 
ence of radioactive 56 Ni to heat the ejecta after initial 
post-explosion expansion of the star. The 56 Ni masses 
inferred from the peak optical brightness of SN-GRBs 
have a broad range, with the brightest, e.g., SN1998bw 
and SN2003dh, requiring several 0.1 M Q . On average, 
however, SN-GRBs are not required to produce more 
56 Ni than the local population of SNe Ic (Soderberg 
2006). In fact, as with low luminosity SNe, e.g., the 
"tailless" SnII, some GRB-SNe may produce little or no 
56 Ni (MacFadyen 2003), as recent observations indicate 
for GRB060505 and GRB060614, two relatively nearby 
(z ~ 0.1) long GRBs with no detected supernova com- 
ponent (Fynbo et al. 2006; Delia Valle et al. 2006; Gal- 
Yam et al. 2006). 

In models of (non-GRB producing) core collapse su- 
pernovae, 56 Ni is produced hydrodynamically in mate- 
rial heated to T > Tm ~ 5 x 10 9 K by the explosion 
shock launched in the core of the star. The amount of 
56 Ni produced depends on the mass inside of the expand- 
ing shock when its temperature declines below Tni- This 
occurs when its radius has expanded to 

where E = E51 x 10 51 erg is the explosion energy and a 
is the radiation constant. The mass inside this radius de- 
pends on the density structure of the progenitor star and 
on how much expansion or contraction occurs before the 
shock reaches a given mass element. In particular, lit- 
tle or no 56 Ni is produced by a shock, even if extremely 
powerful, if it is launched into a low density environ- 
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ment. This may occur if a weak initial explosion ex- 
pands the stellar core so that little mass remains within 
a few 10 8 cm when the strong shock arrives. Produc- 
tion of <~ 0.1M Q of 56 Ni occurs for many pre-supernova 
stars if <~ 10 51 ergs is deposited isotropically by a (quasi- 
) spherical shock on a timescale of ~ Is so that little pre- 
cxpansion of the star occurs before the shock arrives. 
Some of the brightest supernovae, e.g., SN1998bw, re- 
quire energies of up to <~ 10 52 ergs to make the ~ 0.5M Q 
inferred from light-curve modeling. 

The requirement of fast (< Is), isotropic deposition of 
energy for hydrodynamical production of 56 Ni presents 
a serious challenge for models of the SN-GRB central 
engine. First, because the GRB engine must typically 
last 10 s or more for relativistic ejecta to escape the star 
and, second, because GRBs are believed to be highly 
asymmetric explosions. The high degree of beaming and 
long timescale for energy deposition renders collapsar jets 
themselves incapable of producing anywhere near the re- 
quired 56 Ni masses (MacFadyen & Woosley 1999), since 
little mass (< 0.001 M Q ) can be heated to sufficiently 
high temperatures. Therefore, in the original collapsar 
model, with a black hole accretion disk as the central en- 
gine, the 56 Ni is produced in a non-relativistic bi-conical 
wind blown from the disk and constituting a distinct ex- 
plosion component (MacFadyen & Woosley 1999; Mac- 
Fadyen 2003). 

A fundamental problem for the magnetar model, if it is 
to produce a GRB and a supernovae, is the requirement 
that it produce both an isotropic explosion for the 56 Ni 
production and beamed relativistic ejecta. In our model, 
56 Ni can be produced behind a roughly spherical hydro- 
dynamical shock driven by the initial quasi-isotropic ex- 
pansion of the magnetosphere. The expansion becomes 
collimated and the tower formation begins only after the 
stress of the magnetosphere becomes sufficient to bal- 
ance the post-shock pressure. The collimation process of 
the magnetar wind thus involves both a quick isotropic 
expansion followed by a beamed component. We feel 
that this modification to the magnetar scenario, i.e., the 
inclusion of the magnetosphere interaction with the ex- 
terior star, strengthens its viability as a model for the 
long GRB central engine. 

4.3. Restarting the Engine 

We note that the same magnetar can power explosions 
with the degree of collimation that depends on the mag- 
nitude of the outer bounding pressure. A quasi-spherical 
supernova or a highly beamed jet may result from the 
same star at different times as the bounding pressure 
changes. In the collapse and explosion of a massive star, 
the pressure of stellar gas bounding the central magnetar 
may have a complex time history. The star may initially 
collimate the embedded magnetar power into a tightly 
collimated tower responsible for GRB emission. Subse- 
quently, after the star expands and the pressure bound- 
ing the magnetar decreases, the magnetar power will no 
longer be strongly beamed and a normal quasi-spherical 
magnetar outflow will result. Later, however, if material 
not ejected to infinity falls back and accretes, the mag- 
netar will again be surrounded by a bounding pressure 
and its power will be recollimated. X-ray flares observed 
following some GRBs (Burrows et al. 2005; Falcone et al. 
2006; Romano et al. 2006) could result from this process 



(see also, e.g., Proga & Zhang 2006; Perna, Armitage & 
Zhang 2006). 

4.4. Pre-shaping the cavity 

In the previous sections, we have shown that toroidal 
field makes the expanding plasma self-collimating due 
to hoop stress, and propose a spherical cavity with con- 
stant wall properties (i.e. no dependence on polar angle) 
as the simplest model problem. However, the cavity is 
expected in many cases to have lower density near the 
polar axis at fixed radius due to various processes acting 
as the star collapses. Among these are rotational flatten- 
ing and the asymmetric stress from an early magnetized 
wind. First, in order to produce a millisecond magnetar, 
the progenitor star must have been rapidly rotating. We 
therefore expect the collapsed core to be strongly modi- 
fied by rotational effects. In particular the material near 
the rotation axis experiences no centrifugal barrier in- 
hibiting its accretion, resulting in a relatively low density 
in the polar region. A separate effect is that a weaker, 
non-relativistic MHD jet may have been launched along 
the axis earlier, during the collapse of the stellar core 
(LeBlanc & Wilson 1970; Wheeler et al. 2000). In ad- 
dition, an initial MHD wind from the proto-magnctar 
may be concentrated to the poles as in Bucciantini et al. 
(2006). This will push out the cavity in the polar region. 
The subsequent relativistic wind will then expand into 
a cavity pre-shaped by the previous MHD wind. At a 
fixed radius, the pressure and density of the wall will be 
decreased at the poles relative to the equatorial values. 
If these effects arc extreme, the cavity is significantly 
weakened in the polar direction, and a model problem 
consisting of a "magnetar- in-a-tube" is of interest. 

4.5. Pulsar Kicks 

Note that in our picture, most of the magnetically- 
extracted rotational energy of the neutron star travels 
vertically through the two oppositely directed channels. 
Correspondingly, a significant amount of linear momen- 
tum is transported up and down from the neutron star 
and, correspondingly, a back-reaction force is exerted on 
the neutron star from both the top and the bottom. 
The two back-reaction forces are oppositely-directed and 
nearly cancel each other. However, even a slight im- 
balance in the force may have important consequences 
for the overall momentum impacted to the neutron star 
and hence for its terminal velocity. For example, tak- 
ing the total initial rotational energy of the PNS to be 
E ro t — 5 • 10 52 erg, the momentum transported out in 
each direction is P = E Iot /2c ~ 10 42 cgs. Therefore, 
just a 10% imbalance would result in the terminal veloc- 
ity of the neutron star of order of utcrm — 0.1P/Mns ~ 
300 km/sec. 

4.6. Prospects for Numerical Simulations 

In order to gain a solid physical understanding of the 
fundamental physical processes controlling the interac- 
tion of a magnetar with its birth environment, we sug- 
gest a sequence of numerical investigations employing a 
range of well-tested plasma descriptions. Of particular 
usefulness are limiting cases which allow for simplified 
analysis making the key physics more transparent. Sim- 
ulations should cover regions of parameter space where 
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limiting cases overlap with more complete plasma de- 
scriptions. For example, force-free (degenerate) electro- 
dynamics (FFDE) is a useful tool for studying highly 
magnetized plasma for which pressure and inertia are 
negligibly small. In this extreme case, the cavity wall 
would have to be represented by a rigid perfectly con- 
ducting outer boundary condition. While this case may 
not be of general relevance for the realistic physical en- 
vironment, some basic aspects of a bounded rotating 
magnetosphere may be understood using this descrip- 
tion. In addition, the FFDE description has the ad- 
vantage of requiring fewer parameters to specify the ini- 
tial and boundary conditions for a given model problem. 
The results of the time integration can then be more 
easily understood with a minimum of complicating fac- 
tors. Time-dependent force-free codes have already been 
used successfully in the recent years to study pulsar mag- 
netospheres (e.g., Komissarov 2006; McKinney 2006b; 
Spitkovsky 2006). 

It is possible that the full magnetar-in-a-star problem 
can be successfully investigated by a hybrid simulation 
that would employ a relativistic force- free code inside the 
cavity and a relativistic hydrodynamic simulation out- 
side (e.g., R. D. Blandford 2005, private communication, 
McKinney 2006a). 

The next step would be to treat the plasma in the fully 
relativistic MHD regime. There are several relativistic 
MHD codes in existence that have reached the required 
level of maturity (Koide, Shibata & Kudoh 1999; Gam- 
mie, McKinney & Toth 2003; Del Zanna, Bucciantini 
& Londrillo 2003; De Villiers, Hawley & Krolik 2003; 
Fragile 2005; Komissarov 2005; Nishikawa et al. 2005). 
Of interest would be a set of simulations with a range of 
plasma (3. The low-/3 simulations should match on to the 
FFDE case, at least qualitatively. Once these simulations 
are analyzed and the basic physical processes elucidated, 
(3 can be gradually increased enabling an understanding 
of how plasma inertia and pressure affect the dynamics 
of the magnetosphere expansion and collimation. 

The basic process of tower formation and collima- 
tion can initially be explored with two-dimensional ax- 
isymmetric simulations. However, to investigate tower 
stability to non-axisymmetric disruptions, fully three- 
dimensional simulations are necessary. 

Finally, note that the general processes we describe 
here are of interest for many astrophysical systems in- 
cluding non-relativistic central objects (e.g., planetary 
nebulae, see Blackman et al. 2001, Matt et al. 2001). 
For this reason, non-relativistic MHD simulations of 
this problem are of interest in themselves, as well as 
a first step toward fully relativistic MHD. Recent non- 
relativistic MHD simulations indicate that the magnetic 
tower mechanism can operate successfully in a variety of 
astrophysical environments (e.g., Romanova et al. 2004; 
Kato et al. 2004ab; Nakamura et al. 2006, 2007). 

5. CONCLUSIONS 

In this paper we have investigated the millisecond- 
magnetar scenario for the central engine of Gamma-Ray 
Bursts and core-collapse Supernovae. We have focused 
on the interaction between the rapidly-rotating magne- 
tar magnetosphere and the surrounding infalling stellar 
envelope. We have argued that the stellar material pro- 
vides a confining (ram) pressure that has a strong effect 



on both the size and the shape of the magnetosphere. In 
particular, it can channel the highly-magnetized outflow 
originating from the proto-neutron star into two colli- 
mated magnetic towers. 

More specifically, we suggest that the stalled bounce 
shock — a common feature in models of core-collapse 
supernovae — effectively plays a role of a cavity that 
confines the magnetosphere. The cavity's radius, deter- 
mined by the balance between the pressure of the hot 
neutrino-heated gas and the ram pressure of the infalling 
material, stays quasi-stationary at Rq — 200 km during 
the first few hundreds of milliseconds after the bounce. 
To get a feeling for what happens to the magnetar mag- 
netosphere during this stage, we introduce a simplified 
fundamental-physics problem that we call the Pulsar-in- 
a- Cavity problem. A large part of our paper (§ [2]) is 
devoted to investigating this problem. We show that 
since the radius of the cavity is larger than the pulsar 
light-cylinder radius, the magnetic field inside the cav- 
ity continuously winds up. Correspondingly, both the 
toroidal field strength and the magnetic spin-down lumi- 
nosity of the pulsar increase roughly linearly with time. 
The magnetic energy in the cavity then grows quadrat- 
ically with time. We then demonstrate that in the con- 
text of a millisecond magnetar inside a collapsing star the 
magnetic field becomes dynamically important a fraction 
of a second after the bounce. This leads to a subse- 
quent revival of the stalled shock and may result in a 
successful magnetically-driven explosion. As long as the 
expansion of the cavity is non-relativistic, the toroidal 
magnetic field inside it remains larger than the poloidal 
magnetic and electric fields. As a result, the hoop-stress 
collimates the Poynting-flux-dominated outflow into two 
vertical channels that are similar to Lynden-Bell's (1996) 
magnetic towers (see Uzdensky & MacFadyen 2006). 

Finally we discuss the implications of model for several 
observationally-motivated questions relevant to GRBs 
and core-collapse supernovae, such as 56 Ni production, 
late-time X-ray flares, and pulsar kicks. We also outline 
a set of numerical studies that we feel need to be done. 
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